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SUMMARY 


The streamwise flow in a circular cross-section S-duct has been measured 
at four stations for Reynolds numbers of 790 and 48 000 (corresponding to 
Dean numbers of 299 and 18 100) in a water tunnel. The S-duct was formed from 
two 22.5° bends with 48 mm hydraulic diameter and 336 mm mean radius of curva- 
ture. The boundary layer at the inlet to the duct was approximately 25% of the 
hydraulic diameter in the laminar flow: in the turbulent flow it ranged from 
10% to 20% around the periphery of the pipe. Laser-Doppler velocimetry was 
used to measure the streamwise velocity and, in turbulent flow, the associated 
fluctuation. For turbulent flow, the component of velocity at right angles 
to the duct centre-line and the wall pressures are also presented. 

In the first bend, pressure-driven secondary flows are formed but are 
inferred to be larger in the laminar flow case. This is attributed to the 
thicker boundary layer thickness at the inlet. On passing into the second 
bend, the imposed radial pressure gradients give rise to forces which oppose 
the original streamwise vorticity, except for one region near the outer wall 
where the existing streamwise vorticity is sustained, according to the Squire- 
Winter formula for the generation of secondary flow, because of the local sign 
of the radial vorticity. In the turbulent flow the secondary flow on exit from 
the duct and remote from the pipe wall is found to be in the opposite direc- 
tion to that in the first half. 

The core flow, which is defined as fluid with a streamwise velocity great- 
er than 0.90 of the maximum at any given streamwise location, migrates towards 
the outside wall of the first half of the bend and lies towards the inside 
wall of the second half for both Reynolds numbers. Comparison with measure- 
ments in a square cross-sectioned S-duct shows that the core's behaviour is ge- 
nerally similar, with small quantitative differences apparent only in the tur- 
bulent flow. This is because of the smaller secondary flows in the circular 
cross-section duct, attributable to the thinner inlet boundary layers. 

An important objective of this work is to provide observations in suffi- 
cient detail and accuracy for the evaluation of numerical calculation methods. 
Accordingly, the measurements are presented in tabular form and are also avai- 
lable on magnetic tape from NASA Lewis Research Center for ease of reference. 
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1 . INTRODUCTION 


This report presents measurements of the streamwise and one cross-stream 
component of velocity in a S-duct of circular cross-section; it is the sequel 
to an investigation (reference 1) in a similar duct of square cross-section. 
Previous experimental investigations of the flow in S-ducts include observa- 
tions for both fully-developed (reference 2) and thin (reference 3) inlet boun- 
dary layers, mainly in the form of total pressure measurements. In this work, 
the flow at the inlet to the S-duct had thin inlet boundary layers, with the 
purpose of the investigation being the measurement of the velocity characte- 
ristics by a laser-Doppl er velocimeter. An important objective of this report 
is the recording of the velocity field in detail for the evaluation of nume- 
rical iiolutions of the equations of motion. The duct was of relatively mild 
curvature and small centre-line displacement so that there is little likeli- 
hood of flow separation at the walls of the duct. For this reason the flow 
should be particularly amenable to economical 'marching' solution techniques. 
Both laminar and turbulent flow have been investigated because the former 
permits the assessment of the numerical technique alone: the latter tests both 
the numerical technique and any 'turbulence model' which may be incorporated. 

The velocities were measured using a laser-Doppler velocimeter but the 
effect of refraction at the curved duct walls requires a substantially diffe- 
rent experimental procedure from that of reference 1. This is described in 
the following section and is followed by a presentation of the results for 
laminar and turbulent flow. The emphasis of the discussion in section 4 is 
in comparing and contrasting this flow with the results in the square cross- 
sectioned duct. The report closes with a list of concluding remarks. 


2 0 FLOW CONFIGURATION AND EXPERIMENTAL METHOD 


2 0 1 Flow configuration 

The duct, shown schematically in figure 1, was assembled from two 22.5° 
sectors, each of which was manufactured by boring a solid Perspex block 
(n.b. Perspex is a polymethylmethacrylate plastic similar to Plexiglass), The 
ratio of the mean radius of curvature (336 mm) to the hydraulic diameter 
(48,0 + 0.1 mm) was 7.0. The co-ordinate system, also shown in figure 1, is 
identical to that of reference 1. Streamwise location is denoted either by 
distance (X^) along the duct centre-line (expressed in hydraulic diameters, 
with origin at the upstream face of the duct) or by the total angle turned 
through by the duct (9), It is noted that the radial direction always refers 
to the normal to the radius of curvature of the duct centre-line and not to 
the radius of the pipe. 
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The duct was inserted in the water tunnel shown in figure 2, downstream 
of a straight duct of 4.4 hydraulic diameters in length. The water tunnel 
was similar to that used in reference 1, although no boundary layer trip was 
incorporated downstream of the contraction (7.84:1). This was so as to re- 
tain the smallest possible boundary layers on entry to the duct. The experi- 
ments were conducted at Reynolds numbers of 790 and 48 000 (Dean numbers Pf 
299 and 18 100) correspondi ng to bulk velocities (mass flow rate divided oy 
cross-sectional area) of 16.4 rnn s" 1 and 1.00 m s‘l . The mass flow rates we- 
re measured by precision bore flow meters which were insensitive to changes 
in the viscosity of water for the range of temperatures (12°C to 17 0) over 
which experiments were made. 


2 o 2 Laser-Doppler velocimeter, experimental method and measurement tolera nces 

The experimental procedure used in reference 1 for the measurement of 
velocity components lying in the plane of curvature of the bend (observation 
of the velocity at a point in three different directions) could not be appli 
because of the effects of refraction at the curved Perspex/water interface. 
Indeed, cross-stream components of velocity can only be measured along a 
diameter of the pipe: as a result, velocity cross-correlations could not be 
measured. This is not an important restriction^ however, because it is known 
(reference 1) that cross-correlations such as uv are uninfluential except in 
reqions very close to the walls. Surveys of the streamwise component of velo- 
city were therefore made by traversing the optical arrangement shown in figure 
3(a) and applying a ray tracing analysis (given as appendix A) to correct for 
'the effects of refraction on the position of measurement and on the veloci- 
meter's transfer constant (Doppler frequency to velocity). 


Profiles of radial velocity were measured in the direction normal to the 
plane of curvature of the S-duct with the optical arrangement depicted in fi- 
gure 3(b), As for the streamwise measurements, a ray tracing analysis (given 
as Appendix B) was used to correct for the effects of refraction. 


Since the radial component of velocity is always small in comparison 
with the local streamwise component, it is necessary to use ‘frequency shift- 
inq‘ (see reference 4) to generate a measurable Doppler signal. This shift 
was provided by the rotation of a radial diffraction grating (Institute for 
Applied Physics, TN0-TH {Netherlands}, type H) at a nominal speed of 448 
r,p,m,, corresponding to a frequency shift of 0,24 MHz. The relatively large 
value of the shift, as compared to the Doppler frequency corresponding to the 
radial velocities present in the duct, restricted accurate measurement of this 
component to the higher Reynolds number case. 


The optical arrangement of the velocimeter and the electronic signal pro- 
cessing were virtually identical to those used in reference 1 and only the 
principal characteristics of the optical arrangement are summarised in the 
following table, for convenience. 
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Principal optical characteristics of the laser-Doppler velocimeter 


Focal length of 

(a) lens nearest laser 

(b) collimating lens 

(c) imaging lens 

(d) objective lens 

Half angle of intersection (in air) 

Laser light wavelength 
Fringe separation 

Intersection volume dimensions, calculated at 1/e^ 
intensity level (in air) : 

(a) length 

(b) diameter 

Calculated number of fringes in measuring volume 
Photomultiplier pinhole size 


200 mm 
300 mm 
200 mm 
127 mm 

9 o 26 degrees 
632 . 8 nm 
1.97 ym 

1.03 mm 
0 o 17 mm 

85 

0.36 mm 


Transfer constant 

Nominal frequency shift (V-component only) 


0.508 MHz/ (m/s) 
0.24 MHz 


The flows were seeded (see reference 4) with minute dilutions of milk (about 
20 ml in 10^ litres) to increase the scattering-particle concentration. 

Wall static pressures were measured, through tappings let into the wall 
of the duct, by a differential micromanometer. 

Because the experimental procedure in this report is somewhat different 
to that of reference 1, the estimates of the tolerances in the measured quan- 
tities are also different. The following table summarizes these estimates. 
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Ci 


Measurement errors 


Quantity 

r*, z* 

e 


V b (laminar) 

V b (turbulent) 

U 

V 

v 


Systematic error 


Random error 




1 mm 


+ 

0.02 

mm 



0.5 mm 


+ 

0.02 

mm 



1 .5 degrees 


+ 

0.4° 




0.05 

degrees 



nil 




1 % 



+ 

1 % 




nil 



t 

3/4 

% 

up 

to 

0.02 

v b 


+ 

1 % 


up 

to 

0.04 

v b 


+ 

1 % 


up 

to 

0.02 

v b 

up 

to 

3 % 


up 

to 

0.02 

\ 

up 

to 

3 % 


less than 

1 % 



+ 

0.01 

1 


n.b. (a) Quantity '9' refers here to the direction of the local tangent to 
the duct centre-line. 

(b) There is no systematic error in V b (turbulent) because the flow 
meter has been independently calibrated. 

(c) Percentage errors refer to errors in local value. 


The possible errors due to the refraction at the curved surfaces are of 
two kinds. The first is due to inaccuracies in the values of the dimensions 
of the test section, or the refractive indices of the media, which have been 
used in the calculations in appendices A and B. These are likely to be very 
small, however, and unimportant in comparison with the tolerances given in 
the table above. The second is due to residual surface scratches on the Per- 
spex duct walls, remaining after the polishing of the plastic, which might 
reduce the coherence of the laser beams. Since good Doppler signals were 
obtained throughout the experiment, it appears that this source is also negli- 
gi bl e. 

In applying frequency shifting for the measurement of the V-component of 
velocity there are two potential sources of error. The first is associated 
with the accuracy with which the mean rotational speed of the grating _ (which 
is directly proportional to the frequency shift) is known. The rotational 
speed of the grating was measured by frequency counting and is accurate to 
much better than 0.01 V b . The second source is the short term ( r - m - s :' 
stability of the rotation of the grating: this has been estimated as 0.00b v b 
(reference 1) and hence both these sources are negligible compared to those 
already quoted in the table above. 5 



3. RESULTS 


3.1 Laminar flow results 

Flow visualisation (by a combination of the hydrogen bubble technique 
and a plane of laser light) confirmed the laminar, steady nature of the flow. 
Careful observation of the outside wall of the first half of the bend re- 
vealed no evidence of any flow separation due to the adverse pressure gra len 
which exists there. Sample measurements (in both the laminar and turbulent 
flow regimes) on either side of the symmetry plane showed that the flow was 
symmetrical to within the precision of the measurements. 

The streamwise velocities, normalised by the bulk velocity V 5 , are 
shown in contour form in figure 4. There are four stations (1,3,4 and 5) 
which correspond closely with the stations in reference 1. Measurements 
were not made for the location of station 2 in reference 1 since the change, 
relative to station 1, is quite small. The contours have been drawn from 
points linearly interpolated from the measurement grid. In those regions were 
the position of the contour is conjectural, the contour is drawn as a dotted 
line. 


In common with the results of reference 1, the boundary layer thickness 
(defined at 0.95 of the maximum velocity) is about one quarter of the hydrau- 
lic diameter at the beginning of the S-duct. A useful method for comparing 
the results of this work with those of reference 1 is to compare the location 
of the high-speed core of the flow, defined as that fluid which has stream- 
wise velocity greater than, say, 0.90 of the maximum at any given station. 
This core behaves in the same way as for the square cross-sectional bend, 
namely that it is displaced to the outside wall in the first half of the b- 
duct due to the induced secondary motions and remains at r - 0.3 in the 
second half of the duct. By the outlet of the bend there is a large simu- 
lation of low speed fluid near the outside wall of the second bend, which 
was also noted in reference 1, with an extensive region (r ^ 0.5, z ~ • ) 

where aU/az* is positive. As noted in reference 1 and 2 the importance of 
this region in the second half is that it provides a source of secondary 
(i.e. streamwise vorticity) motion which is of the same sense as that esta- 
blished by the first half of the S-duct and thus continues the accumulation 
of low speed fluid which was begun in the first halfof the S-duct. ihis 
is in contrast to regions where aU/az* is negative (in other words, the wa 
boundary layers) which generate streamwise vorticity of a sense which atte- 
nuates the secondary flow set up by the first half of the S-duct. 


The measurements are tabulated in appendix C. 

3.2 Turbulent flow results 

Flow visualisation of the turbulent flow was also carried out at the out 
side wall of the first half of the S-duct. Again, no evidence of flow sepa- 
ration was observed although it should be stated that observation of this re- 
gion is substantially more difficult because of the high flow velocities. 
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The streamwise isotachs for the turbulent flow are depicted in figure 5, 
with values normalised by the bulk velocity, V b . The boundary layers are, 
as expected, thinner than for the laminar flow and vary. from 0.10 to 0.20 
of the hydraulic diameter around the periphery of the pipe. This range is 
slightly thinner than that of the turbulent flow in the square cross-sec- 
tioned S-duct of reference 1, because of the shorter upstream tangent and 
the absence of a boundary layer trip. The core flow at station 1 is dis- 
placed towards the inside of the first half of the S-duct more markedly 
than for the corresponding station of the square cross-sectioned S-duct 
( re f erence i). Although the core migrates towards the centre of the pipe 
by station 3, there is little evidence of any accumulation of low speed fluid 
at the inside wall, indicating weak secondary motion. The boundary layer in 
the round duct remains slightly thinner than in the square duct at. the same 
station. On passing into the second half of the S-duct, the core is again 
displaced towards the inside surface and this displacement persists to the 
exit plane of the S-duct. The persistence of this displacement is, once 
again, more pronounced than in the case of the square cross-sectioned duct 
(reference 1). As in the case of the laminar flow, a region of positive 
aU/3z* exists near the outside of the second half of the duct. which strength- 
ens the existing secondary motion in accumulating low streamwise. velocity 
fluid. It can be seen, for example, by noting the relative positions of the 
U = 0.7 V b contour, that low-speed fluid is more extensively distributed in 
the case of the circular, rather than square, cross-sectioned duct. 

The profiles of V-c omponent velocity in the z* direction (for r* = 0.5) 
are given in figure 6. The magnitudes at the entry to the duct are extre- 
mely small (about 0.02 V b ) and in the direction of the displacement of the 
core of the streamwise isotachs. By the inflection plane (station 3) compo- 
nents of about -0.05 V b are observed in the core flow region with a measured 
maximum of about + 0.12 V b near the wall. These values are smaller than 
those (-0.09 V b and + 0.13 V b ) in the square-sectioned duct. In the second 
half of the duct the magnitudes decrease, and in the core flow eventually 
reverse, as a result of the change in the direction of the curvature of the 
S-duct. The importance of this change is well illustrated by the centre-line 
profile of radial velocity (figure 7), which clearly shows the abrupt alte- 
ration in the development on passing from the first to the second half of 
the duct. In reference 1 it was observed that negative components of V-velo- 
city occur near the wall at station 5, but this is not evident here. This 
is consistent with the generally smaller secondary flows in this duct, al- 
though it is possible that negative components exist at locations closer to 
the wall than the extent of the measurements. It should also be recalled 
that the magnitudes of velocity which are under consideration are comparable 
to the systematic error of the velocimeter. 

The streamwise variation in the wall static pressure around the pipe is 
presented in figure 8. The form of the graphs is similar to, if somewhat 
smaller than, those occurring in the square-sectioned S-duct of reference 1. 

The turbulence levels for the streamwise component of velocity are shown 
as contour levels in figure 9 and as cross-stream (z*) profiles for the v 
component in figure 6. The turbulence level in the core is low 
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(u ~ v - 0.02 V^), as would be expected of the unstrained part of the flow, 
and higher values occur within the wall boundary layers. Qualitative agree- 
ment with reference 1 is good, although the 'pocket* of comparatively high 
turbulence which is observed at r* = 0.8, z* * 0.4 for stations 4 and 5 is 
absent from the square cross-sectioned duct. Its occurrence here is the 
corollary of the more extensive region of low momentum fluid in the vicinity 
of the outside of the second bend (comp, above). 

The mean and r.m.s. velocities are tabulated in appendix D as table 1 
and the wall pressure measurements as table II. 


4. DISCUSSION 

Developing flows in ducts can be discussed, at least qualitatively, in 
terms of streamwise acceleration and deceleration, brought about by the im- 
pressed streamwise pressure distribution, and cross-stream convection of 
streamwise momentum, due to secondary flows. These two mechanisms are, 
in many cases, competing ones in that the pressure distribution acts so as 
to accelerate the flow on the inside of a bend whereas streamwise vorticity 
(i.e. secondary flow) convects the low-speed boundary layer to the inside of 
the bend, thereby displacing the high speed core away from this region. The 
exclusion of the effects of shearing stresses (whether laminar or turbulent) 
from the discussion is not to deny their importance (see, for example, the 
numerical experiments of reference 6) but their influence does not introduce 
any qualitatively new features. 

A simple analytical expression for the generation of streamwise vorti- 
city is (reference 7) 



<5o) s = 2 0)^60 

(1) 

where 

_ 3W 3V 
w s 3r*‘ 3Z 

(2a) 

and 

_ 3W 3U 

W r 30 3Z 

(2b) 


are the streamwise and radial components of vorticity, respectively, W 
is the velocity component in the z direction and r* is the radial co-ordinate 
directed towards the centre of curvature in each bend. Since 3W/30 « 3U/3z, 
equation (1) states that the magnitude of the streamwise vorticity which is 
generated depends on the value of 3U/3z or, what is equivalent in this case, 
the boundary layer thickness. In the laminar flow, this thickness was similar 
for both circular and square cross-sectioned S-ducts so that the behaviour of 
the core flow is virtually identical. In contrast, the boundary layer in the 
turbulent flow was thinner at the inlet to the circular cross-sectioned S-duct, 
with the result that the generation of streamwise vorticity is also smaller. 
Hence the flow is more influenced by the streamwise pressure gradient than by 
secondary flows and thus the behaviour of the core flow resembles that of a 
potential flow more closely than was the case for the square cross-sectioned 
S-duct. 
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5. CONCLUDING REMARKS 


1 This report presents, in detail, the velocity characteristics of laminar 
and turbulent developing flow in a S-duct formed from two 22.5° sectors of 
a bend of circular cross-section, with a radius ratio of 7.0. 

2. Pressure-dri ven secondary flows arise in the first half of the S-duct with 
a measured maximum of about 0.12 of the bulk velocity in turbulent flow. 

On passing through the second half of the S-duct the secondary flow esta- 
blished upstream is attenuated near the pipe walls, and reversed in the 
core flow region. 

3. In common with observation in the square cross-sectioned S-duct, there is 
a region near the outside wall of the second half of the S-duct where the 
local sign of the radial vorticity results in a secondary flow which sus^ 
tains that established upstream. This leads to a substantial local redis- 
tribution of the streamwise isotachs. 

4 The laminar flow development is similar in both the circular and square 
cross-sectioned ducts. In the turbulent flow, the core flow in the case 
of the circular cross-sectioned S-duct lies closer to the inside wall of 
each half of the S-duct. It is argued that this is the result of thinner 
boundary layers at the entrance to the duct, as compared with the square 
cross-sectioned duct. 

5. This report presents benchmark measurements against which the results of 
numerical calculation methods can be compared. 


Note: The author's names are in alphabetical order. 
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APPENDIX A 


DETERMINATION OF BEAM PATH FOR STREAMWISE 
VELOCITY MEASUREMENTS 


A 1. PROBLEM STATEMENT 

The position of the intersection of the two laser beams, and the angle 
with which the intersection occurs, are functions of the traverse co-ordina- 
tes because of the effects of refraction at the Perspex/air and Perspex/water 
interfaces. The half-angle of intersection (B) is a quantity of interest 
because of the dependence of the transfer constant (K) of the velocimeter on 
this angle: 

U(t) = Kf (t) A 1(a) 

where K E liTHW A1(b) 

Here U(t) and f(t) are the instantaneous values of streamwise velocity compo- 
nent and Doppler frequency, respectively and X and n are the wavelength of 
laser light and the refractive index of the medium in which the beam inter- 
section takes place. The purpose of this Appendix is to determine, by emu- 
lation, the position of the intersection and the angle 3 as a function of the 
traverse co-ordinates of the optical bench. The analysis is derived below in 
detail for the first half of the S-duct and the changes which are required 
for the second half are briefly described under a separate heading. 


A 2. ANALYSIS FOR THE FIRST HALF OF THE S-DUCT 

Figure A1 shows the attitude of the laser-Doppler velocimeter (LDV) for 
the measurement of streamwise velocity at any streamwise station within the 
first half of the S-duct. Let the (orthogonal) traverse directions of the op- 
tical bench be x and y, with origin as shown; both of these quantities are 
measured. The origin of co-ordinates for position relative to the duct, r 
and z, is also shown. In plan view, the optical axis of the LDV strikes the 
duct symmetrically and thus it is sufficient to determine the path of one beam 
only. The beams enter the Perspex wall at point A and are refracted at this 
interface to point B, the Perspex/ water interface. Beam intersection takes 
place within the water at point C. The analysis proceeds by considering the 
problem in three stages: refraction at the air/Perspex interface, refraction 
at the Perspex/water interface and finally the beam path within water. 

Refraction at the air/Perspex interface at point A 

The beam is incident at angle a, found as follows. Let F be the virtual 
focus of the unrefracted laser beam and 6 the point on the axis of the gene- 
rator of the S-duct which contains point A. Considering triangle AGF, the 
sine rule gives 
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A 2 


sin(a) _ sin(7T-(f>) 


where R a is the radius of the outer air/Perspex interface and $ is the half- 
angle of beam intersection in air. Hence angle a can be found and the angle 
of refraction, b, can be found from Snell's law 

n a sin(a) = npSin(b) A 3 

where n a and np are the refractive indices of air and Perspex respectively. 
Refraction at the Perspex/ water interface at point B 

The refraction at this point is best analysed using vector algebra. The 
components of these vectors are resolved relative to a dextral system of co- 
ordinates, shown in figures A 1 and A 2 , with origin at point B with i_, j_ and 
J< being the unit vectors in the direction of the axes of the co-ordinate sys- 
tem. 

The direction of the incident beam is given by p (a unit vector). Its 
direction is determined by the angle Y. This angle Ts found from its comple- 
mentary angle, u, by application of the sine rule to triangle ABG: 



Y = TT-U 


A 4(a) 


sin(u) 

R a ' 

sin(b) 

BG 

A 4(b) 

with the length BG, 

from figure A 1 , 

and A 3 



BG = R + 

(£) cos(e) 

A 4(c) 

where 

e = arc 

sin (^ 2 *), from figure A3 

A 4 (d) 


Here R is the radius of curvature of the centre-line of the S-duct and d is the 
diameter of the pipe. 

Hence £ = cos(J -Y)i_ + cos(y) j_ + 0 k A 5 


Refraction takes place about the normal to the Perspex/water interface at 
point B, with direction given by r (a unit vector) increasing from the incident 
to the transmitting medium. Its direction is defined by angle e (equation A 4 
(d)) and hence 


r_ - 0 £ + cos(e)£ - cos(J -£ )j< A 6 

The direction of the refracted beam is denoted by £ (a unit vector) and is 
set by refraction: it is a required quantity. Let xn denote the angle of inci- 
12 



dence: its value is given by 


by definition of the 
unit vectors. If x w 
law. 


cos(x p ) = £ * I 

vector scalar product and noting that both p and r_ are 
denotes the angle of refraction then, applying Snell s 

" w sin( v = v in( y At 


which allows X w to be found. The direction of 4 can now I be found.! since the 
laws of refraction require that the refracted ray lies in the plane of inci 
dence. The vectorial representation of this statement is (reference , 
p. 65, equation 4.8) 


n..q = (n,..cos(x w ) - ncosfx ))r + n 6 


A 9 


Beam path within water 

It is convenient for the subsequent analysis to find the direction of 

the line BD (figure A 1). Let this be denoted by s a All bv the 
in the plane defined by i and j., and is determine? (figures A 2, A!) by the 

angle BDG (angle n). From figure A1 

A 10(a) 

n = p- t ' ' 

^ , . , A 10(b) 

where p + a - 4> 

(external angle of triangle AFG equals the sum of interior opposites) 

and S + b = T A10(C) 

(external angle of triangle ABG equals the sum of interior opposites). Since 
angles a, b and Y are known from the preceding, angle n is also known. 

Hence 


s_ = cos(n)i + cos(-?j- - n)i + 0 k 


All 


The half-angle of intersection, which is one of the required quantities of 
the analysis, can now be found. This angle, BCD, is denoted as e: referring to 
figure A 4, if the angle DBC is a, then 


3 = 


- a 


But a is the angle between vectors q and s_ and is given by 

cos (a) = £ • s_ 


A 12 


A 13 


and hence the angle $ can be found. 
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The other quantity which is required is the position of the point of in- 
tersection (point C) relative to the origin of co-ordinates r and z, figure 
Al. 

r co-ordinate . From figure Al, 

r = (R + |) - (DG - CD cos(a)) A 14(a) 

where a is the angle between line CD and the i_-j_ plane and is as yet unknown. 
Now 

DG = BG cos(n) A 14(b) 

from figure Al, with BG and n given by equations A(4c) and A 10(a) above. 

The length CD is 

CD = BD tan (a) A 14(c) 

from figure A 4, with a known from equation A 13 and 

BD = BG sin(n) A 14(d) 

The angle a is found as follows. From Figure A 4, the direction of line CD 
can be written as 

t = _q - s_ cos(a) A 15(a ) 

The vector defining the normal to the i_- j_ plane is k_ and, from the definition 
of angle a 

k • t = | kj |t|cos(| -a) A 15(b) 

where |kjsl (unit vector) and |tj can be found from A 15(a). Hence o can be 
found. 

z co-ordinate . From figure A 1 

z = x + CD sin(a) A 1 6 

Since x, CD (equation A 14(c)) and a (equation A 15(b)) are all known, z can 
be found. This completes the analysis for the first half of the duct. 
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A3. ANALYSIS FOR THE SECOND HALF OF THE S-DUCT 


The analysis is exactly the same, 
half of the duct. Accordingly, only a 
the changes. 


mutatis mutandis, as for the first 
brief statement is made concerning 


Refraction at the air/Perspex interface at point A 

Referring to figure A5, equation A2 for the angle a becomes 

sin(ir-a) _ sin (4> ) 

'~*T 


A 17 


where Ri is the radius of the inner air/Perspex interface. The equation de 
fining angle b, equation A3, is unaltered. 

Refraction at the Perspex/ water interface at point B 

Equations A4(a), (b) and (c) defining angles y, u and the length BG 
are re-written as follows: 


sinhO _ sin(iT-b) 
R, BG 


BG 


= R - (j) cos ( e ) 


A 18(a) 
A 18(b) 


For this analysis, the angle u is unnecessary and e is given by equation A 4(d), 
as before. The refraction at point B is represented by figure A 6 & nd b y 
inspection the representation of vectors £ and r, equations A 5 and A 6, are 
unaffected, as are the equations for X p and xw Tequations A 7 and A 8). Hence 
vector q, defined by equation A 9, stands. 


reads 


Beam path within water 

Equation A 10(a), for angle n, is unaltered, but the equation for p now 

a = P + <j> A19 (*) 

from triangle AFG and the equation for c now is 

b = ? + Y A19 ( b ) 

from triangle ABG. The direction of vector £ is 

£ = cos (n) i_ - cos(tj- - n)j_ + Ok. A 20 

and the change in sign in the j component, relative to equation Ml, 1S 
and is evident in figure A 6. ^he equations for ^(equation A 12) and a (equa 
tion A 13) remain unaltered for the second half of the duct. 



r co-ordinate . From figure A 5 


r = (DG + CD • cos (c )) - (R - |) A 21 

which replaces equation A 14(a). The equations for lengths DG, CD and BD 
(A 14(b), (c) and (d) ) are unaffected, as is the analysis for finding the 
angle a (equations A15(a) and (b)). 

z co-ordinate . Equation A 16 holds for this case. 


A 4 . TYPICAL RESULTS 

For the determination of the beam path in the S-duct, the following va- 
lues were used: 


Refractive index of air, n a 1.00 

Refractive index of water, n w 1.33 

Refractive index of Perspex, n p 1.490 

Angle of intersection in air, <f> 9.26 degrees 

Radius of curvature of outer surface, R a 377.70 mm 

Radius of curvature of inner surface, R-j 293.97 mm 

Radius of curvature of bend centre-line, R 336.00 mm 

Diameter of pipe, d 48.00 mm 


The two graphs in figure A7 depict the change in the transfer constant 
(relative to that in air) of the LDV, aK, and the position of the beam inter- 
section, r, both as a function of the radial traverse co-ordinate, y, for the 
first half of the duct. Both graphs compare the values for zero elevation 
(z* = 0.0) with those obtaining at z* = -0.63 and show that the influence of 
elevation on AK(y) and r(y) is extremely small. 
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APPENDIX B 


DETERMINATION OF BEAM PATH FOR r-COMPONENT 
VELOCITY MEASUREMENTS 


B 1. PROBLEM STATEMENT 

The position of the intersection of the two laser beams of the laser- 
Doppler veloci meter (LDV), and the angle with which this intersection takes 
place, are functions of the traverse co-ordinate in the same way as for ap- 
pendix A. It is the purpose of this appendix to determine these two quanti- 
ties by calculation. The problem is restricted to the case where the axis 
of symmetry of the LDV coincides with the vertical diameter of the pipe: 
that is the r co-ordinate is always r=d/2, where d is the diameter of the 
pipe. The analysis which follows is based on that of reference 

B 2. ANALYSIS 

Fiqure B1 shows the attitude of the LDV for the measurement of r-compo- 
nent velocity at any streamwise location of the S-duct. The traverse distan- 
ce of the optical bench is t, with origin of traverse taken as that position 
at which the beam intersection occurs at the Perspex/water interface. The 
optical axis of the LDV strikes the cross-section symmetrically and it is the- 
refore sufficient to determine the path of one beam only. The beams enter the 
Perspex wall at point A (or A') and are refracted to point B (or B), the Per- 
spex/water interface: beam intersection takes place at point c within the wa 

ter. Analysis proceeds by considering the displacement of point E (the vir- 
tual focus of beam A' B’ ) given a traverse movement t and from this calculates 
the position of point C. 

Refraction at the air/ Perspex interface 

Let D be the virtual focus of the undisturbed laser beam when the traverse 
is positioned such that intersection occurs at the Perspex/ water interface. 
Suppose the optical bench to be raised by a distance t as shown in thefigure. 
The beam path is now A' B' and E' is the virtual focus of the beam as it pro- 
pagates through the Perspex: it is required to find the distance BE , denoted 
as x. If the angle FBG is a right angle, then 


and 


FB = AA l 
A'A= t*tan(<{>) 


B 1 


where <$> is the half-angle of beam intersection in air and which is known from 
measurement. The angle FE'B is equal to the angle (a) of refraction at point 
A' which, by Snell's law, is given by 


n a sin(4> ) = n p sin(a) 


B 2 


where n_ 

From triangle 


a and n D are the refractive indeces of air and Perspex respectively. 
FE'B 


FB _ t . tanU ) 
tanfa) tan (a) 


B 3 
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Refraction at the Perspex/ water interface 


Two cases are distinguished. If 


B 4 


where d is the diameter of the pipe then the beam path is as shown in figure 
Bl. The beam is incident at an angle b to the normal and is found by apply- 
ing the sine rule to B'E'G: 


that is 


E'G = B'G 
sin(b) sin (it - a) 

( ?~ x ) _ (d/2) 
sin(b) sin(Tr-a) 


B 5(a) 
B 5(b) 


The angle of refraction is c and is given by Snell's law 


n p * sin(b) = n w * sin(c) 


B 6 


The half-angle of intersection of the beams, 3, which is one of the required 
quantities, is 


3 = angle (B'E'B) + angle (C B ' E ' ) 


B 7 (a) 


since external angle of triangle B'CE' is equal to the sum of interior opposi- 
tes. Since 


and 

then 


angle (B'E'B) = a 

angle (CB'E* ) = c-b 
3 = a+c-b 


B 7 (b) 
B 7(c) 
B 7(d) 


The location (z) of the point (C) of the intersection of the beams can now be 
found by application of the sine rule to triangle B'CG. 


Ul _ (d/2) 
sin(c) sin(Tr-3) 

where |z| denotes the magnitude of z. 


B 8 


If 


x > £ 


B 9 


then the beam path is as shown in figure B2. Equation B 5(a), from which the 
angle b is found, is now 


(x-d/2 ) _ (d/2) 
sin(b) si n (a ) 


BIO 


The angle of refraction, c, is unaltered (equation B6). The half angle of 
18 



intersection is 


3 = angle (GE'B) - angle (E'B'C) 
angle (GE'B) = a 
angle (E'B'C) = c-b 
hence 3 = a-(c-b) 

The location of point C is now given by 

z = (d/2 ) B 1 2 

sin(c) sin 3 


Bll(a) 
B 11 (b) 
B 11 (c) 
B 11(d) 


B 3. TYPICAL RESULTS 

The graph in figure B3 depicts the change in the transfer constant (re- 
lative to that in air for <j> = 9.11 degrees) of the LDV, AK, as a function of 
the position of the measuring volume, z*. The change in this quantity is. 
much greater than was the case for the streamwise measurements, detailed in 
Appendix A and is partly due to the larger angles of incidence which occur. 
However it is noted that even when the measuring volume lies at the centre 
of the duct (z* = 0), and hence there is no refraction at the Perspex/water 
interface at point B' , there is an increase of about 12% in K over the value 
obtaining in air. This is due to the change (from 632.8 nm in vacuo) in wave- 
length of the laser light because of the refractive index of water (n w - 1.33) 
as well as refraction at the air/Perspex interface at A'. 
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APPENDIX C 


TABULATED DATA 

22.5°/22. 5° S-DUCT, CIRCULAR CROSS-SECTION 
LAMINAR FLOW STREAMWISE VELOCITIES (REYNOLDS 
NUMBER = 790, DEAN NUMBER = 299, BULK VELOCITY 
= 16.4 m/s) 


Notes: (1) 'THETA' is angle of turning of bend, 
measured from S-duct inlet. 

(2) ' R 1 and 'Z' are cross-stream co- 
ordinate directions, defined in 
figure 1 . 

(3) 'U ' is the streamwise component of 
vel ocity. 
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APPENDIX D 


TABULATED DATA 

22. 5°/22.5° S-DUCT, CIRCULAR CROSS-SECTION 

TURBULENT FLOW : REYNOLDS NUMBER = 48 OOO 
DEAN NUMBER = 18 100 

BULK VELOCITY = 1.00 m/s 

TABLE I 

MEAN & r.m.s. COMPONENTS OF STREAMWISE AND 
RADIAL VELOCITIES 

Notes: (1) 'THETA' is angle of turning of bend, 
measured from S-duct inlet. 

(2) * R * and 'Z ' are cross-stream co-ordinate 
directions, defined in figure 1. 

(3) ' U ‘ is the time-mean streamwise component 
of velocity and 'UD ' is the corresponding 
r.m.s. velocity fluctuation 

(4) 'V' is the time-mean radial component 

of velocity and ' VD ' is the corresponding 
r.m.s. velocity fluctuation 

(5) 'XH' is distance along duct's centre-line 
(in hydraulic diameters) with origin at 
the duct inlet. 
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TABLE II 

WALL PRESSURE DISTRIBUTION 
Note : 'Pressure coefficient' is defined as 

C p . ( p - p ref> 

p = — r~2 — 
i pV b 

_3 

where p is water density (1000 kg m ) 
is bulk velocity (1.00 m/s) 


STPEA^WI St 
STATION 
UH> 

RADIAL 
POSIT ION 

m 

SPANsmt 

POS I II ON 
1 1 > 

PRtr sy»t 
c jefficient 

- 1.00 

0.0 

0 .0 

0.0 (PRFFJ 


0.5 

-1 .0 



1.0 

0.0 

- 0.015 

1.10 

0.0 

0.0 

0.070 


0 . £■ 

-1 .0 

- LJ .0 50 


1.0 

0.0 

- 0.203 

2.20 

0,0 

0.0 

3.0 A 2 


o . 5 

-1 .0 

- 0.070 

2.38 

1.0 

0,0 

- 0,203 

3.30 

0 . 0 

U .0 

- 0 . 27 A 


0.5 

“ 1 .0 

- 0.137 


1.0 

C .0 

- 0.038 

A . AO 

0.0 

0.0 

- 0.326 


0.5 

- 1.0 

- 0.153 


1.0 

0.0 

- 0.053 

6.50 

0,0 

0,0 

- 0.196 


0.5 

- 1.0 

- U . 20 A 


1.0 

0.0 

- 0.190 
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APPENDIX E 


Roman 

S 

De 


DEFINITION OF SYMBOLS 

Characters 

Pressure coefficient (figure 8) 
Hydraulic diameter (48 mm) 


ref 


R 

Re 

U 

0 

V 

v b 

v 

W 

X H 

Z 

z* 


Dean number: De = 


ID 


X 


. Re 


Pressure at wall 

Reference value of P (X^ = 0.0, r* = 0, z* = 0) 

Radial co-ordinate direction (figure 1) 

Radius of curvature of suction surface 
Radius of curvature of pressure surface 
Normalised radial co-ordinate: r* = ^ r o 

V r o 

Mean radius of curvature: R = i ( r -j + r 0 ) 

Reynolds number based on and D^ 

A, 

Velocity in streamwise direction 

R.m.s. fluctuating velocity in streamwise direction 

Velocity in radial direction 

Bulk mean velocity (mass flow rate/ duct area) 

R.m.s. fluctuating velocity in radial (r) direction 
Velocity in spanwise direction 

Distance along duct centre-line, expressed in hydraulic 
diameters (origin lies at 0 = 0°) 

Spanwise co-ordinate direction (figure 1) 

Duct ha If -width (24 mm) 

Normalised spanwise co-ordinate: z* = z__ (see figure 1) 

z. 



Greek Characters 

0 Angle of turning in streamwise direction (see figure 1) 

<P Ha If -angle of beam intersection 

w $ Streamwise component of mean vorticity (equation 1,2a) 

Radial component of mean vorticity (equation 1,2b) 
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Figure 2 Plan view of water tunnel 
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Figure 3(b) Optical arrangement of LDV for radial velocity measurements 








Station 4, X H =3 85, &= 22 5° + 9° 


Station 5 , X H = 5.22, 0- = 22 5° + 20 25 


Figure 5 Turbulent flow: isotachs of U/V b at successive streamwise stations 
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Figure 6 Turbulent flow: profiles of V/V^ and v/V^ for r * = 0.5 (n.b.V is positive in the positive 
r direction) 
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Figure 8 Turbulent flow: wall pressure distribut 


Station 4 , X H = 3.85, 6- = 22 5°+ 9.0° 


Station 5, X H = 5 .22 , 9- = 22 .5°* 20 .25° 


Figure 9 Turbulent flow: contours of u/V^ (x 100) at successive streamwise stations 
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Figure A 2 Refraction of beam at point B for the first half of S-duct. 

Refer also to the definition of the vector origin, figure A1 
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Figure A3 Section across duct, passing through point B 
(both first and second halves of S-duct) 
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Figure A7b. - Calculated radial displacement (r) of beam intersection 
as a function of traverse co-ordinate, y. 
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OPTICAL AXIS 
OF L.BV, 


OPTICAL axis 
OF L..D.V 



Figure B 1 Path of laser beam for V-component 
velocity measurement. Intersection 
volume formed below plane of symmetry 
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Figure B2 Path of laser beam for V-component 
velocity measurement. Intersection 
volume formed above plane of symmetry. 
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